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Abstract 

In this article, we consider the compressible Navier-Stokes equation 
with density dependent viscosity coefficients. We focus on the case 
where those coefficients vanish on vacuum. We prove the stability of 
weak solutions for periodic domain fl — T N as well as the whole space 
il = M. N , when N = 2 and N = 3. The pressure is given by p = p 1 , 
and our result holds for any 7 > 1. In particular, we prove the stability 
of weak solutions of the Saint- Venant model for shallow water. 

1 Introduction 

This paper is devoted to the Cauchy problem of the compressible Navier- 
Stokes equation with viscosity coefficients vanishing on vacuum. Let p(t, x) 
and u(t, x) denote the density and the velocity of an isentropic compressible 
viscous fluid (as usual, p is a non-negative function and u is a vector valued 
function, both defined on a subset f2 of Mr). Then, the Navier-Stokes 
equation for isentropic compressible viscous fluids reads (see |LL59j ): 

d t p + div (pu) = . . 

d t (pu) + div (pu®u) +V x p - div (hD(u)) -V(gdivu) = ^ ' 

where p(p) = p 1 , 7 > 1, denotes the pressure, D(u) = l[Vu + *Vu] is the 
strain tensor and h and g are the two Lame viscosity coefficients (depending 
on the density p) satisfying 

h>0 h + Ng>0 (2) 
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(h is sometime called the shear viscosity of the fluid, while g is usually re- 
ferred to as the second viscosity coefficient). One of the major difficulty of 
compressible fluid mechanics is to deal with vacuum. The problem of exis- 
tence of global solution in time for Navier-Stokes equations was addressed in 
one dimension for smooth enough data by Kazhikov and Shelukhin KS77 , 
and for discontinuous one, but still with densities away from zero, by Serre 
Scr8Q] and Hoff |Hof87j . Those results have been generalized to higher 
dimensions by Matsumura and Nishida |MN79j for smooth data close to 
equilibrium and by Hoff Hof 95b| . |Hof 95a in the case of discontinuous data. 

Concerning large initial data, Lions showed in lLio98l the global existence 
of weak solutions for 7 > 3/2 for N = 2 and 7 > 9/5 for N = 3. This result 
has been extended later by Feireisl, Novotny, and Petzeltova to the range 
7 > 3/2 in FNP01 , and very recently by Feireisl to the full system of 
the Navier-Stokes equations involving the energy equation lFei04j . Other 
results provide the full range 7 > 1 under symmetries assumptions on the 
initial datum (see for instance Jiang and Zhang £ JZ03 ). All those results do 
not require to be far from the vacuum. However they rely strongly on the 
assumption that the viscosity coefficients are bounded below by a positive 
constant. This non physical assumption allows to get some estimates on the 
gradient of the velocity field. 

The main difficulty when dealing with vanishing viscosity coefficients 
on vacuum is that the velocity cannot even be defined when the density 
vanishes. The first result handling this difficulty is due to Bresch, Desjardins 
and Lin |BDLQ3j . They showed the L 1 stability of weak solutions for the 
following Korteweg's system of equations: 

d t p + div (pu) = . . 

dt(pu) + div (pu (£> u) + V x p — z^div (pD(u)) = upVAp. 

The result was later improved by Bresch and Desjardins in BD03 to include 
the case of vanishing capillarity (k = 0), but with an additional quadratic 
friction term rp\u\u (see also BD02 ). The key point in those papers is to 
show that the structure of the diffusion term provides some regularity for 
the density thanks to a new entropy inequality. However, those estimates 
are not enough to treat the case without capillarity and friction effects k = 
and r = (which corresponds to equation (^Q) with h(p) = p and g(p) = 0). 

The main difficulty, to prove the stability of the solutions of is to 
pass to the limit in the term p(u®u) (which requires the strong convergence 
of ^fpu). Note that this is easy when the viscosity coefficients are bounded 
below by a positive constant. On the other hand, the new bounds on the 
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gradient of the density make the control of the pressure term far simpler 
than in the case of constant viscosity coefficients. 

Our result is in the same spirit as the one of Bresch, Desjardins and Lin 
and makes use of the same entropy inequality, first discovered by Bresch 
and Desjardins in BD02 for the particular case where h(p) = p and g{p) = 
0. We actually use a slightly more general estimate, which holds for any 
viscosity coefficients h(p),g(p) satisfying the relation: 

g{p) = ph'{p)-h{p). (4) 

This estimate first appeared in a Note by Bresch and Desjardins B D04| in 
the context of Korteweg systems of equations. However, we will see that 
the capillary term is by no means necessary to the derivation of the crucial 
estimates which thus hold for the compressible Navier-Stokes system 

Our main contribution is to show the L 1 stability of weak solutions of 
(PQ) under some conditions on the viscosity coefficients (including @) but 
without any additional regularizing terms. The interest of our result lie 
primarily in the fact that our conditions allow for viscosity coefficients that 
vanish on the vacuum set. It includes the case h(p) = p, g(p) = (when 
N = 2 and 7 = 2, we recover the Saint Venant model for Shallow water), but 
our conditions on h and g will exclude the case of constant viscosity h(p) = p, 
d(p) = £■ Indeed, it is readily seen that @ implies that g{p) = £ = — p, and 
thus p + ^ = 0. In this border line case we thus lose all informations on the 
derivatives of u. It is worth pointing out that while we can gain regularity 
on the density with this new estimate, we have to loose regularity on the 
velocity (on the vacuum set). 

Note that the main difficulty will be to establish the compactness of yfpu 
in L 2 strong, and the key ingredient to achieve this is an additional estimate 
which bounds ^fpu in L°°(0, T; L 2+2a (ft)) for some small a > (the usual 
entropy estimate only gives a bound in L°°(0, T; L 2 (ft))). 

For the sake of simplicity we will consider the case ft = IR^ and the case 
of bounded domain with periodic boundary conditions, namely ft = T N . For 
the same reason we consider only power pressure laws although the result 
could be extend to non monotonic pressure law of the form of |Feif)2| . Note 
that the result holds for any power 7 > 1 under appropriate assumptions on 
h and g. Classically, L 1 stability is considered as the main step to prove the 
existence of weak solutions. To obtain the existence of weak solutions, one is 
thus left with the technical task of constructing a sequence of approximated 
solutions verifying the a priori estimates. Although this final step is in 
most cases quite standard, we point out that in this particular situation it 
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seems highly non trivial because of the complexity of the additional entropy 
inequality 

In the next section, we state the assumptions on the viscosity coefficients, 
define precisely the notion of "weak solutions" and state our main results. 
In Section |3J we recall the well known physical energy inequality and state 
the key estimates. The proof of Theorem 12.11 is detailed in Section 21 For 
the sake of completeness, we give in Section [S] the proof of the entropy 
inequality of Bresch and Desjardins in the context of compressible Navier- 
Stokes equation. 



2 Notations and main result 

Let Q denote a subset of Mr. We assume that Q is either the whole space 
M N or a bounded domain with periodic boundary conditions (O = T N ). 
For the sake of simplicity, we will take D(u) = Vu, though the full strain 
tensor could be considered without any additional difficulty. This leads to 
the following system of equations: 

dtp + div (pu) = (5) 
d t (pu) +div(>u® u) + Vxp 1 - div (h(p) Vu) - V(g(p)divu) = 0,(6) 

with initial conditions 

p\t=o = p o >0, pu\t=o = m a . (7) 

Before introducing the notion of weak solution, let us state the assumptions 
we make on the viscosity coefficients. 



Conditions on h(p) and g(p): 
First we consider f(p),g(p) verifying: 

g{p)=ph'{p)-h{p). (8) 

As stated in the introduction, this structure constraint is fundamental to 
get more regularity on the density. Moreover, we assume that there exists a 
positive constant v G (0, 1) such that 

ti(p)>v, h(0)>0 (9) 

\g'(p)\ < \h'{ P ) (10) 

vh(p)<h(p)+Ng(p)<-h(p). (11) 



4 



When 7 > 3 and N = 3, we also require that 

liminf-^->0, (12) 
for some small e > 0. Let us make some remarks about those assumptions. 
Remark 2.1 The functions 

Hp) = P, g{p) = o 

satisfy In fact, any linear combination of p k with k > 1 is an 

admissible function for h(p). 

Remark 2.2 The lower estimate in Ml]) is trivial when g > 0, while the 
upper estimate is trivial when g < 0. Together they yield: 

\g{p)\ < C u h(p) V P >0. 

This inequality and il(J\) will be necessary to pass to the limit in the term 
V{g{p n )divu n ). 

Remark 2.3 Condition makes the proof simpler, but is not optimal. 
However, condition Ml\) is necessary to control the viscosity term and to- 
gether with (0), it yields 

N-l + u ^ h\p) N-l + l/v , „ 
-^p—^Wj- N P ' forall P >0, 

Cp (N-l)/N+u/N < h{p) < Cp (N-l)/N + l/(Nu) p > : 

Cp (N-l)/N+l/(Nu) < h ^ < Cp (N-l)/N+u/N t p < 1 ( U ) 

In particular, we must have h(0) = 0. Moreover, this shows that if we do not 
assume the "best" h(p) we can take is h(p) = pi N ~ l )/ N + u / N . This is 
actually enough to prove the stability of weak solutions for all 7 when N = 2 
and for 7 < 3/2 when N = 3. However, if we assume h(p) ~ Cp 2 / 3+u for 
small p and h(p) ~ Cp for large p, then we can take any 7 G (1,3) when 
N = 3. 



and so 
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Notion of weak solutions 

We say that (p,u) is a weak solution of Q5I6JI on 17 x [0,T], with initial 
conditions (JJJ) if 

p g L°°(o,T,L l (n) n^(o)), 

^uGL°°(0,T; (L 2 ^))^), 

MP)V« G L\0,T; (w£>\Sl)) N * N ), <?(p)div U G L 2 (0, T; W^ 1 ' 1 ^)), 
with p > and (p, *Jpu) satisfying 

f d t p + div {yfpyfpu) = . n v , 

\ p{0,x) = p a (x) 

and if the following equality holds for all (p(t,x) smooth test function with 
compact support such that <p(T, •) = 0: 

/ m • 99(0, •) dx + / / *Jp(yfpu)d t (p + y/pu <g) y'pii : Vy? 
Join 

rp 

+ [[ p 1 div ipdx - (h(p)Vu, V^> - (g(p) (div u) , (div if)) = 0,(14) 
Join 

where the diffusion terms make sense when written as 
(h(p)Vu, Vip) = 

Hp), 



{y/puj)duipj dx dt — I (^/p~Uj)2h' (p)di^/pdiipj dx dt, 



(g(p) (div u) , (div ip)) 
9(P), 



and 



_ (yfp~Ui)dijipj dxdt-J (yfp~Ui)2g (p)d iy /pdjipj dxdt. 

In particular, the fact that the diffusion term h(p)\7u (and g(p)divu) lies in 
L 2 (0,T; (W 1 ~ 1,1 (J2)) nxn ) will follow from the fact that 

h'(p)V^GL°°(0,T;Ll c (n)), and h(p)/ */p G L°°(0, T; L 2 oc (J7)), 

and similar conditions on fir(p). This will be provided by assumptions l|10[). 
© and (H3I). 



Main result: 

The main result of this paper is the following: 
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Theorem 2.1 Assume that 7 > 1 and that h(p) and g(p) are two C 2 func- 
tions of p satisfying conditions H&\)- \11\) (together with \12\) if 7 > 3 and 
N = 3). Let (p 

n> u n ) n £fq be a sequence of weak solutions of J5H6J) satisfying 
entropy inequalities U8\) . \21\) and h25\) . with initial data 

Pn\t=0 = Po( x ) and PnU n \t=0 = = Po ( x ) u o ( x ) > 

where p™ and are such that 

P™>0, p n ^Po inL\n), pX-> PoUoinL 1 ^), (15) 
and satisfy the following bounds (with C constant independent on n): 

/* I 77, 1 2 1 /* 1 

/ Po^r + TPVdx < C, / -\Vh(p n )\ 2 dx < C, (16) 

Jn 2 7 - 1 A} Po 

and 

f |„.n|2+(5 

p n o l -^—dx<C, (17) 

/or some small 5 > 0. 

Then, up to a subsequence, (p n , yfp~^u n ) converges strongly to a weak so- 
lution of (Ep-© satisfying entropy inequalities ilfy). \21\) and \25\) (the den- 
sity p n converges strongly in C°((0, T); L^{Q)), ^fp^u n converges strongly 
in L 2 (0, T; L 2 oc (£l)) and the momentum m n = p n u n converges strongly in 
L\Q : T-L\ oc ^l)), for anyT>0). 



3 Entropy inequalities and a priori estimates 

In this section, we recall the well-known energy inequality and state the 
main inequalities that we will use throughout the proof of Theorem 12.11 

The usual energy inequality associated with the system of equations (|B1l6|) 
can be written as: 



d f u 2 1 I' f 

— / p— + p^dx+ / h(p)\Vu\ 2 dx+ / g(p)(divu) 2 dx<0. (18) 

dtj 2 7-1 J J 



This inequality can be established for smooth solutions of (l">l()|) by multi- 
plying the momentum equation by u. 

When h and g satisfies h(p) + Ng(p) > and if the initial data are taken 
in such a way that 



f u 2 1 

£0= Po-^r H tpI dx < +00, 

J n 2 7-I 
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then (HBJ) yields: 



\WP~ u \\L°°(0,T;L 2 (n)) < Cj 
l|p||L°°(0,T;L7(n)) < C 



(19) 



Furthermore, Hypothesis (jlljl gives: 



1 1 \/M^) Vu| |£,2(o,T;£,2(n)) < C- 



(20) 



Finally, integrating Q with respect to x yields the natural L 1 estimate: 



Unfortunately, it is a well-known fact that those estimates are not enough 
to prove the stability of the solutions of Q5I6JI . In particular, the fact that 
p 1 is bounded in L°°(0, T; L X (S7)) does not implies that pZ converges to p 7 . 

However, further estimates can be obtained by mean of the following 
lemma (the proof of which is postponed to Section EJ: 

Lemma 3.1 Assume that h(p) and g(p) are two C 2 functions such that 
10) holds true. Then, the following inequality holds for smooth solutions of 



with <p such that 

<p' = -. (22) 
P 

This lemma is similar to the result of D. Bresch and B. Desjardin in BD04 , 
in which the same inequality was derived when capillary effects are taken 
into account. 

We immediately see that since the viscosity coefficient h(p) is an increas- 
ing function of p and when the initial data satisfies 



PUL^iOtT;!, 1 ^)) 



< c. 



JH1). 




(21) 




inequality ((2Tj) yields: 
1 

-||^Vt£(/))|| LO o (0iT;i 2 (n)) = \\h (p)V v / p|k°°(o,T;i, 2 (n)) < C> ( 23 ) 



and 



\Wh'{p)p^Vp\\ L 2^ T . L 2 m < C. 



(24) 



S 



Under assumption Q on h, those estimates will give additional control on 
the density p and on the pressure /? 7 , which will be enough to prove the 
stability of weak solution. 



Finally, we shall make use of the following result: 
Lemma 3.2 Assume 

h(p) + Ng(p) > vh(p) 

for some v G (0,1) (which is a part of ill}) )- and let 5 G (0,^/4). Then, 
smooth solutions of J3H2J) satisfy the following inequality: 

d ' M2+5 -dx + ^ J h(p)\u\ S \Vu\ 2 dx 

p\u\ 2 dx I 



dt / 2 + 5 





where \Vu\ 2 = £^ J2j \d%u 



2 

, 3 \ . 



This inequality is quite simple to establish and will be essential in the proof 
of Theorem 12. II to prove that ^/p^u n is bounded in L°°(0, T; L 2+2a (Q)) (see 
Lemma 14 .3JI . Note, however, that to derive further estimates from this 
inequality, we need to control the right hand side of ()25[). Inequality (|T%|) 
immediately provide a bound on J p\u\ 2 dx, so the problem will be to control 

/ 27— 6/2 \ 2/(2— 5) 

enough power of p to get a bound on J ( p ^4 j dx. This will be 

achieved using (|24[). Of course, we also need to assume that the initial 
condition satisfies 

Po—^ — dx < C. 

Proof of Lemma \3.°A Let 5 G (0,^/4). Multiplying (JBJ) by we get: 

\u\ 2+s , /■ „M 2+<5 



dz + pu-V—dx 



+ y /i(p)|ii| 5 (Vii) 2 (ix + 5 j h(p)\u\ s 2 UiUkdjUidjUk dx 
+ / g(p)\u\ s (div u) 2 dx + 5 j g(p)\uf~ 2 UkUjdiUidjUk dx 



+ I \u\ s u ■ Vp 1 dx = 0. 
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Since 

(divu) 2 = d i u i d 3 u 3 < \( diV % + 9 i u ^ - N \ Vu ^ 

i 3 i 3 

condition (|TT|) yields: 

|2+<5 r U,|2+5 



u 



2 + 5 



dx + / pu-V 



2 + 5 



dx + v h(p)\u\ 6 (Vu) 2 dx 



+ j \u\ 5 u-Vp~<dx 
— fij h(p)\u\ S ~ 2 UiUkdjUidjUk dx 

+5 g(p)\u\ s ~ 2 u k UjdiUidjU k dx, 



and since 5 < z//4, we deduce: 

1 2+5 



\u\ + f \u\ + V f 

pdt^j dx + J pu- V^-y dx + - J h(p)\u\ s (Vu) 2 dx 

u\ 6 u ■ Vp^dx < 0. 



Moreover, multiplying (|5j) by 1 2 ' +<? and integrating by parts, we have 



|2+5 



2 + 5 



1 2+5 



dtp dx — / pu - V— dx = 



2 + 5 



and summing the last two inequalities, we get: 



/I 1 2+5 f f 

P-j^j dx + \\ Hp)H S \Vu\ 2 dx< \u\ s u ■ V/> 7 dx 



d 

37 / P 



dt J r 2 + 5 

It remains to bound the right hand side. We have: 
\u\ s u ■ Vp 7 dx 

< (VN + 5) 



p" f \u\ s divu dx — 5 j p^\u\ & 2 u{u ■ X7)udx 



< (VN + 5) ( I h(p)\u\ d \Vu\' z dx 



1/2 



< V - I h(p)\u\ 5 \Vu\ 2 dx + C u I ^\u\ s dx, 
4 J J h(p) 



f-. \u\ s dx 

Hp) 
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where the last term satisfies (if S £ (0, 2)): 




and the lemma follows. 



We now have all the necessary tools to prove Theorem 12. II 



4 Proof of Theorem 12.11 

We now present the proof of Theorem 12. 11 To begin with, we need to make 
precise the assumptions on the initial data. 

Initial data: 

We recall that the initial data must satisfy (|16|) . and 1|17|) to make use of all 
the inequalities presented in the previous section: 

p n Q is bounded in L 1 n L 7 (ft), p n Q > a.e. in Q 
p n \u n \ 2 = \m n \ 2 /p n is bounded in L 1 ^) 

Jo~2V<p(p2) = Vh{p n )/Jp~2 is bounded in L 2 (fl), (26) 



r \ u n \ 2 + s 
/ p2 — ^— — dx < C for some small 8. 



With those assumptions, and using inequalities (fT8]> and (|2Tj) . we deduce the 
following estimates, which we shall use throughout the proof of Theorem 12. II 

\WPn U n\\L°°(0,T);L 2 {n)) < C 

l|P n|U°°(Q ,T;L 1 nZ,i , (n)) — & (27) 
1 1 V KPn) V«n| |i2(0,T;i 2 (n)) < C 
and _ 

ll^ / (Pn)Vy / P n | |L°°(0,T;L 2 (n)) < C 
ll\/^ / (Pn)pn~ 2 V/J rl || L 2( 0)T;L 2 (C )) < C 



In view of hypothesis on the viscosity coefficient Q, the bounds (|27|) 
and (ESJ) yields: 

||v / Pn^ n nlU 2 (0 ) T;i 2 (Q)) < C 

l|VVp n |U-(o,r ; La(n)) <C (29) 

l|Vp^ /2 ||L2(0,T;L 2 (n)) < C 
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The proof of Theorem 12.11 will be divided in 6 steps. In the first two 
steps, we show the convergence of the density and the pressure (note that 
the convergence of the pressure is straighforward here) . The key argument of 
the proof is presented in the third step: We prove that ^/p~^u n is bounded in 
a space better than L°°(0, T; L 2 (Q)). In turn, this will give the convergence 
of the momentum (step 4) and finally the strong convergence of ^/p~^u n in 
Lf oc ((0,T) x O) (step 5). The last step adresses the convergence of the 
diffusion terms; It is mainly technical and of minor interest. 

Step 1: Convergence of ^fp~^. 

Lemma 4.1 If 'h satisfies J^, then 

y/p^i is bounded in L°°(0, T; H 1 ^)) 
dt\/7>n is bounded in L 2 (0, T; if _1 (f2)). 

As a consequence, up to a subsequence, -^/pn converges almost everywhere 
and strongly in L 2 (0,T; Lf oc (Q)). We write 

\f~fhi — ► ^fp a.e and Lf oc ((0,T) x $7) strong. 
Moreover, p n converges to p in C°(0, T; L^(fi)). 

Proof. The second estimate in (|29|) . together with the conservation of mass 
\\Pn(t)\\ L i(n) = WPn^WL^n) gives the L°°(0, T; H 1 ^)) bound. Next, we 
notice that 

dty/p^, = -^y/p^divu n - u n ■ V^/p^ 
= ^ \/p^div u n - div (u n ^fp^) 

which yields the second estimate and, thanks to Aubin's Lemma, gives the 
strong convergence in L^ oc ((0, T) x Q). 

Sobolev imbedding insures that yj~p~^ is bounded in L°°(0, T; L 9 (0)) for 
q G [2, +oo[ if N = 2 and q E [2, 6] if N = 3. In either cases we deduce that 
p n is bounded in L°°(0, T; L 3 (J7)), and therefore 



Pri^n — V Pny Priori 



is bounded in L°° (0, T; L 3/2 (Q) ) . 
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The continuity equation thus yields dtp n bounded in L°°(0, T, W 1,3 / 2 (f2)). 
Moreover, since Vp„, = 2y / p^Vy / p^', we also have that Vp n is bounded in 

L°°(0,T;L 3 / 2 (ft)), hence the compactness of p n in C([0,T];L^(Q)). 



Step 2: Convergence of the pressure 

Lemma 4.2 The pressure pZ is bounded in L 5 / 3 ((0, T) x VL) when N = 3 
and i r ((0, T) x fl) for all r G [1, 2[ w/ien N = 2. In particular, pZ converges 
to p 1 strongly in Lj oc ((0,T) x f2). 

Proo/. Inequalities (ESI and (EZj) yield G L 2 (0,T; H 1 ^)). 

When JV = 2, we deduce G L 2 (0, T; for all q G [2, oo[. So 

is bounded in L x (0, T; LP(fi)) n L°°(0, T; L 1 ^)) for all p G [l,oo[, hence ^ 
is bounded in L r '((0,T) x 0) for all r G [1,2[. 

When N = 3, we only get G L 2 (0, T; L 6 (Q)), or 

Since /9n is bounded in L°°(0, T; L 1 (Q)), Holder inequality gives 

IITII I I T I I I I T I 1^' ^ <f 

\\Pn\\L 5 / 3 ((0,T)xn) - H^nllLoo( 0! T;Li(n)) N^nllH(o,T;L3(n)) - °" 

hence ^ is bounded in L 5//3 ((0, T) x fi). 

Since we already know that pZ converges almost everywhere to p 7 , those 
bounds yield the strong convergence of pZ in Lj oc ((0, T) x Q). 



Step 3: Bounds for ^fp~^u n 

Lemma 4.3 If ^ < 3, or if N = 3, 7 > 3 and holds, then 

yfp^u n is bounded in L°°(0, T; L 2+2a (S7)) 
for some small a > 



This Lemma is really the corner stone of the stability result. As a matter 
of fact, at this point, the main difficulty is to prove the strong convergence 
of sfp^Un in L 1 (0, T; L 2 oc (r2)). A first consequence of Lemma 14.31 is that it 
will be enough to prove the convergence almost everywhere. However, since 
we are only able to prove the convergence of the momentum p n u n (see Step 
4, which makes use of Lemma 14.31 as well), we need to control yfp~^u n on 
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the vacuum set {p(t, x) = 0} (and prove that it converges to zero almost 
everywhere). And this fact also will be a consequence of Lemma |4.3I (see 
Step 5). 

Proof. The proof of Lemma 14.31 relies on Lemma 13.21 for 5 small enough 
(5 £ (0,z//4)), we have: 



d f u , v f, , lSl „ l2 , 
— / p— — -dx + - / hip) u d Vd dx 
dtj 1 2 + 5 2 J yHn 1 1 1 

(2-0/2 



< I / I cix ( / pH a da? 1 .(30) 



h{p) 

Using l|27jl. we deduce: 

2/(2-5) \ ( 2 ~ <5 )/ 2 



|u |2+5 / r / p 2 7 -<5/2 



Condition @ yields h(p) > i^p and so 



(/(^— ) 2/M -) M/2 



Using Lemma 14.21 we readily check that the right hand side is bounded L 1 
in time (for small 5), without any condition when N = 2, and when N = 3 
under the condition that 

5 

27 - K g7. 

which gives rise to the restriction 7 < 3. In either cases, we deduce 

dtj H 2 + 5 

and ljl7JI gives the lemma. When N = 3 and 7 > 3 we need the extra 
hypothesis (|12j) to achieve the same result. 



Finally, for a < 5/2, we have: 

h «(/-*-.)" (/,~) ! 
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with q = (1— (2 + 2a) /(2 + <5)) , so that the exponent of p goes to 1 when 
a goes to zero. In particular, it is less that 3 for a small enough, and since 
p n is bounded in L°°(0, T; L 3 (il)), we deduce Lemma POl 

Step 4: Convergence of the momentum 

Lemma 4.4 Up to a subsequence, the momentum m n = p n u n converqes 
strongly in L 2 (0, T; L^ e (Q)) (for some positive e) and almost everywhere 
to some m(x, t). 

Note that we can already define u(x, t) = m(x, t) / p(x, t) outside the vacuum 
set {p(x,t) = 0}, but we do not know yet whether m(x,t) is zero on the 
vacuum set. 
Proof. We have 

Priori — y/ Pn^J Pn^m 

where ^fp^ is bounded in L°°(0, T; L q (Q)) for q G [2,+oo[ if N = 2 and 
q € [2,6] if N = 3; Since ^/p^u n is bounded in L°°(0, T; L 2 (Q)), we deduce 
that 

p n u n is bounded in L°°(0, T, L g (Q)) for all q G [1,3/2]. 
Next, we have 




Using Lemma f4.3l and (|29|) . it is readily seen that the second term is bounded 
in L°°(0, T; L 1+£ (r2)) for some small e > 0, while the first term is bounded 
in L 2 (0,T,L q {n)) for all q G [1,3/2]. Hence 

V{p n u n ) is bounded in L 2 (0, T; L 1+e (Q)). 

In particular, we have 

Pn u n bounded in L 2 (0, T; W 1,1+£ (fl)). 

It remains to show that for every compact set K C fi,we have 

d t {p n u n ) is bounded in L 5 / 3 (0, T; W~ 2 ' 3/2 (K)). (31) 
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As a matter of fact, we observe that Wq' 3 (K) C L 1+1 ' £ (K) for small e (for 
N = 2 or 3), and therefore 

L 1+£ (K) C W- l ^ /2 {K) C W- 2 ' 3/2 (K), 

so (|31|) together with Aubin's Lemma, yields the compactness of /? n ii n in 
i 2 (0,T;L 1+E W). 

To prove (|31|) . we use the momentum equation ijB]). first noticing from 
Lemma 14.21 and Lemma 14, 31 that 

div(^u„®^u„) e L°°(0, T; W -1,1+e (if)) 
Vp^ e L 5 / 3 (0,T; W" -1 ' 1+e (ir)), 

So we only have to check that V(h(p n )Vu n ) and V (g(p n ) div u n ) are 
bounded in L°°(0,T; W~ 2,3 / 2 (K)). To that purpose, we write 



h(p n )Vu n = V(h(p n )u n ) - u n Vh(p n 
(and similarly with g{p n ))- The second term in ()32|) is 

V/i(p„ 



(32) 



U n Vh(p r , 



IPn 



2yfJ^U n h! {p n )V yfp~ r , 



which is bounded in L°°(0, T; L 1+£ (0)) thanks to (J2EJ) and Lemma Ol The 
first term in ()32|) can be rewritten 



h(Pr, 



'Priori 



which is bounded in L°°(0,T;W 1,3 / 2 (S7)) thanks to the following lemma: 

Lemma 4.5 For all compact set K , h(p n )/y/p^ and g(p n ) / \fp~a o.re bounded 
in L°°(0,T;L 6 (K)). 

The proof of this Lemma is a bit technical in full generality and will be 
postponed to AppendixlAl However, note that, in the particular case h(p) = 
up, we have h{p n )/yfp~^ = yfp and Lemma 14.51 follows straightforwardly 
from Lemma 14.11 

We deduce that h(p n )Vu n and g(p n )div u n are bounded in 
L°°(0, T; W' 1 ' 3 / 2 ^) + L 1+£ (K)), 
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and since L l+S (K)) C W~^' 3 ' 2 (K) we can conclude that h(p n )Vu n and 
g(p n )divu n are bounded in L°°(0, T; W~ 1,3 ^ 2 (K)), which conclude the proof 
of Lemma 14.41 



Step 5: Convergence of y/p^u n 

Lemma 4.6 The quantity y/p^u n converges strongly in L 1 and L 2 oc ((0, T) x 
0) to m/y/p (defined to be zero when m = 0). 

In particular, we have m(x,t) = a.e. on {p(x,t) = 0} and there exists 
a function u(x, t) such that 

m(x,t) = p(x,t)u(x,t) 

(note that u is not uniquely defined on the vacuum set {p(x,t) = 0}). 

Proof. First of all, since m n /y/p~^ is bounded in L°°(0, T; L 2 (Q)), Fatou's 
lemma yields 

2 

liminf — — dx < oo. 

Pn 

In particular, we have m(x,t) = a.e. in {p(x,t) = 0}, and if we define 
m 2 /p to be when m = 0, we have 

2 

m 

— dx < oo. 
P 

Lemma IP1 implies that yffh\u n \ is bounded in L°°(0, T; L 2+2a (£l)) for a 
small a. It is thus enough to prove the convergence almost everywhere, or 
in Lj oc ((0,T) x fi), to prove the strong convergence in L 2 oc . 

First of all, we note that in {p(x,t) / 0}, y/p^u n converges almost 
everywhere to m/y/p. So, if we denote the vacuum set by 

A = {p(x,t) = 0}, 

we deduce 

VfaUnlGA > —^ l £A a - e -- ( 33 ) 

\JP 

To controle y/p^u n on the vacuum set, we introduce the set 

Bl, = {Pl! {2+S) \un\ > M}, 
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for M > 0. We then cut the L 1 norm as follows: 



in 



f ... + f / 



Fpu^n — \dxdt = I 

The L°°(0, T; L 2 (£2)) bound and (|33|) gives the convergence of the first 
integral: 

(34) 



hBV, \A\VPn U n — \dxdt 



0. 



Moreover, Lemma 14.31 and Tchebychev's inequality yields 

C 



\t*M\ S M2 : 



and so 



B" 



Priori 



m/^p\<^h\{J 



Pn\Un\ 2 + \m\ 2 /pdx 



1/2 



< 



C_ 
M' 



(35) 



Finally, on ZB^ n A, we have 

Iv^nl < M P y 2 - 1 /(2+<5) _^ o a.e. , 

since p n — > a.e. and 1/2 — 1/(2 + 5) > 0. So lcs« n^lvTv^nl converges 
almost everywhere to 0. In particular, the L°°(0, T; L 2 (f2)) bound yields 



Since we defined m/yfp to be on A, we also have 



in 



hence 



/ 



m 



ics™ nA\\/~Pn u n - ^j=\dxdt 
VP 



(36) 



Putting (jHU, © and © together, we deduce 



lim sup 

n— >oo 



m C 
— \dxdt < 

VP 



M 



for all M > 0, and so ^fp~^u n converges to ^= in L 1 ((0, T) x O) strong. The 
lemma follows. 
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Step 6: Convergence of the diffusion terms 
Lemma 4.7 We have 

h(p n )Vu n — > h(p)Vu in V' 

and 

g{p n )divu n — ► g(p)divu in V' 
Proof. Let <fi be a test function, then 



h(p n )V u n (j) dx dt = — j h(p n )u n V (f) dx dt + j u n Vh(p n )<j) dx dt 
''''' r p^u n V '(j> dx dt + / \fp~^u n — ^p^-V p n (f> dx dt 



i— V r»o lb T — 1 I V r ib" iv i— 

\fP n J 

Thanks to Lemma we know that is bounded in L°°(0, T; Lf oc (Q)). 

Moreover, since h(p n )/y/p < VyJ~p~^, this term converges almost everywhere 
to h(p)/y/p (defined to be zero on the vacuum set). Therefore, it converges 
strongly in L 2 OC ((0, T) x Q); This is enough to prove the convergence of the 
first term. 

Next, we note that 

k ' ipn) -Vp n = W( Pn ) 



y/Pn 

with ip'(p) =h'{p)lyfp= y/p<p'(p). Since 

\V^(p)\ 2 dx = f p\V^(p)\ 2 dx, 



we have that Vijj(p n ) is bounded in L°°(0, T, L 2 (Q)). Moreover, (fT3)l yields 
h'{p) < Cp~ 1/2+u/3 when p < 1 

and so 

tp{p) < Cp u/3 when p < 1. 

Therefore, an argument similar to the proof of Lemma 14.51 shows that if>(p n ) 
is bounded in L°°(0, T; Lf oc (Q,)). Since it converges almost everywhere (tp is 
a continuous function), it converges strongly in L 2 oc ((0, T) x O). It follows 
that 

V^(p n ) ^ V^(p) Lf oc ((0, T) x ft)-weak . 

A similar argument holds for g(p n )div u n using the fact that \g(p)\ < Ch{p) 
and |</(p)| < C7i'(p). 
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5 Proof of Lemma 13.11 

We conclude this paper by giving the proof of the estimate (|21|) . To that 
purpose, we have to evaluate 



d_ 

Tit 



^p\u\ 2 + pu-Vip(p) + ^p\Vip(p)\ 2 



d /* 1 

dx + — -p 1 dx. 

dtj 7-1 



Step 1: First of all, we recall the usual entropy equality: 

dx = — J h(p)\Vu\ 2 dx — J g(p)\divu\ 2 dx 



d_ 

dt 



2 n 1 7- V 



Step 2: Next, © gives 



pdt dx — / div pudx 

Z i z 



pVu : Vip(p) ® Vip(p) dx + I p 2 ip'(p)A(p(p)divudx 
+ J p[V(p(p)] 2 div udx 



and so 



j p N^P)\ dx = _ j pVu . v^(p) ® Vtp(p) dx 
+ / p 2 (p'(p)Aip(p)divudx 



+ J p[Vip(p)] 2 div udx (37) 



Step 3: It remains to evaluate the derivative of the cross-product: 
— j pu-V(p(p)dx = j V V(p) " dt(pu) dx + y pu ■ dtV f(p) dx 

Vtp(p) ■ dt(pu) dx — J div (pu)(p' (p) dtp dx 
V<p(p) ■ d t (pu) dx+ (div (pu)) 2 ip'(p) dx(38) 
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Multiplying © by Vip(p), we get: 
J Vp(p) ■ d t (pu) dx 

(h(p) + g(p))A(p(p)dW udx + J Vn : Vcp(p) <g> Vh(p) dx 
Vip(p) ■ Vh(p)div udx — j Vp{p) • Vp 7 dx 
V<p(p)div (pu (8) u) dx, 
where we used the fact that 

J ^(g(p)divu) ■ Vip(p) dx = — J g(p)Aip(p)divudx 

and 

div (h(p)Vu) ■ Vp(p) dx 

dj(h(p)djUi))diip(p) dx 
di(h(p)djUi))dj(p(p) dx 

dih(p)djUidjip(p) dx — J diUidjh(p)dj<p(p) dx 

diiiih(p)djjtp(p) dx 
Vu : Vh(p) (8) Vp(p) dx — j Vh(p) ■ Vp(p)div u dx 
h(p) Ap(p)d\v u dx 



Step 4: When ip, h and g satisfies © and J22J), then (J37J) and 
yields 

j t ^Jpu ■ Vtp(p) + p^Mt dx} + f V<p<j>) • Vpdx 

V<y?(/?)div (pu (g) u) <ix + / (div (pu)) 2 dx. 
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Finally, we have 

V p(p)div (pu®u)dx + j p'(p)(drv (pu)) 2 dx 

-p'(p)u ■ Vpdiv(pu) — (p'(p)Vp(pu ■ Vu) + (p' (p)(div pu) 2 dx 
ptp'(p) div u div (pu) — p<p'(p)V p(u ■ Vu) dx 
p 2 (p'(p) (div u) 2 + ptp' (p)u ■ Vpdivu — pip' (p)\7 p(u ■ X7u) dx 
so using ((2*2*)) and (JHJ), we get 

V(/?(p) div (pu®u) dx + j p' (p)(di\ (pu)) 2 dx 

ph'(p)(divu) 2 + V(/i(p)) • udivu - V(h(p))(u ■ X7u) dx 
ph'(p)(divu) 2 — h(p)(div u) 2 — h(p)u ■ Vdivudx 
+ J h(p)diUjdjUi + h(p)u ■ Vdivudx 
= J (ph' — h)(divu) 2 + h(p)diUjdjUi dx 
= J g(p) (divu) 2 dx + j h(p)djUidiUj dx 

which yields 

j t U pu ■ Vp(p) + p \^}l dx} + J Vp(p) ■ Vp' dx 
< / g(p)(divu) 2 dx + / h(p)\ Vu\ 2 dx, 



and the proof is complete. 

A Proof of Lemma 14.51 

We shall only prove the result for h(p n ) / \fp^. Using the fact that 
\g(p)\ < Ch{p), and \g'(p)\ < Ch'(p) for all p, 
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a similar proof follows for g{p n ) / \fpn 
Note that In view of Q13|). we have 



Hp) . n v ., . -, 
< Gp if p < 1, 



so we only need to control t= for large p n . This will be achieved differently 



VP 

h(pn) 

depending on the dimension. 

When N = 2, the fact that ^ is bounded in L°°(0,T; H 1 (Q)) and 
Sobolev's inequalities implies that p n is bounded in L°°(0, T; for all 

p € [l,oo [. Moreover, in view of (|13j) . we have 

%) < / Cp 1 ^ if p > 1 



Vp - I Cp" if p< 1 ' 

So there exists q D > 1 such that h(p n )/^/p~^ is bounded in L°°(0, T; L q (Q)) 
for all (/ > q Q . In particular, h{p n )/yfp~^ is bounded in L°°(0, T; L P (K)) for 
all p G [1, 00 [ for any compact set K. 

When N = 3, we note that 
and since conditions (JHJ) and lfTT|) yields 

^V = 9(p ) + Mrt >M£0+M>^ (A 



we have 



iv (hi^y< C \h f (p)v^p\. 



So inequality (|23j) yields 



||V (^j=-^j \ \L°°(0,T;L 2 (n)) < C (39) 

When ft = R 3 , Sobolev's inequalities implies that h(p n )/y/~p^ is bounded in 
L°°(0, T; L 6 (ft)). When ft is a subset of M 3 , we note that JT3J) gives 

Up) r cpV 6 + 3 /^ if p > 1 

"VP " I Cp 1/6+i//3 if p < 1 ' 
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So there exists a constant s < 1 such that 
'h(p n ) 



IPn 



1 G L°°(0, T; L (f2)) 



Moreover 

'h(p n ) 



V 



1 



< 



V 



HPn) 
h(p 



V 



IPn 



G L°°(0,T;L 2 (ft)), 



using the fact that s — 1 < 0. It follows that (h(p n ) / \/Pn) s ^p n >i is bounded 
in L°°(0, T; which in turn gives 



(^VWieL~(0,T;L a (n)), 



for all si G (s,3s). As long as 3s < 1, we can repeat this argument with 3s 
instead of s. Eventually, this will lead to 



h(p r , 



l Pn >i GL°°(0,T;L 2 (O)), 



which, together with (|39j) implies that (h(p n )/ y/p^)l Pn >i is bounded in 
L°°(0,T;L 6 (ft)). 
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